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Abstract. We study the dynamics of a prolate spheroidal microswimmer in Poiseuille flow for different 
flow geometries. When moving between two parallel plates or in a cylindrical microchannel, the swimmer 
performs either periodic swinging or periodic tumbling motion. Although the trajectories of spherical and 
elongated swimmers are qualitatively similar, the swinging and tumbling frequency strongly depends on 
the aspect ratio of the swimmer. In channels with reduced symmetry the swimmers perform quasiperiodic 
motion which we demonstrate explicitely for swimming in a channel with elliptical cross section. 

PACS. 47.63.Gd Swimming microorganisms - 47.15.G- Low- Reynolds- number (creeping) flows - 02. 30. Ik 
Integrable systems 



1 Introduction 



Swimming on the micron scale has attracted a lot of at- 
tention among physicists, since they try to understand 
the fascinating strategies, microorganisms employ to over- 
come the constraints of low Reynolds number hydrody- 
namics [ll[2]. On the other hand, artificial microscopic 
swimmers or active particles with various locomotion mech- 
anisms have been constructed recently also with the goal 
of using controllable environments to study principal prop- 
erties of their diffusive and collective motion, as well as 
locomotion under external fields such as gravity or shear 
flow [3,4,5 ,6 ,7,8 ,9, 10 , 11 . The experiments are accompa- 
nied by detailed theoretical studies analyzing hydrody- 
namic (see e.g. p!2 l[T3l[T4l [T5 l [T6 l [T7 1[T8ll 1 9( 120] ) . stochastic 
(see e.g. [2T | l22 l[23l [2T ) . chemotactic (see e.g. [25ll26]). op- 
tical [2^ or thermophoretic [28 effects on the dynamics 
of active particles. 

In nature, microorganisms move in an aqueous envi- 
ronment and often have to respond to fluid flow like mi- 
croscopic plankton in the sea [29 , pathogens in the blood 
stream ^30^, or sperm cells in the fallopian tubes when 
swimming towards the egg [31]. Rheotaxis of sperm cells 
[32] or bacteria [3T allow directed motion along flow gradi- 
ents. Bottom-heavy microorganisms show stable orienta- 
tions in flow resulting in layering of phytoplankton in the 
ocean ^34j or hydrodynamic focussing in Poiseuille flow 
[351136]. Neutrally buoyant microswimmers, however, fol- 
low periodic trajectories when swimming in Poiseuille flow 
[37 |[38| r3Q j . Finally, microswimmers in time-dependent flow 
easily show chaotic motion and non-trivial transport phe- 
nomena [391140] . 



In our previous Letter we have studied the dynamics 
of spherical microswimmers in cylindrical Poiseuille flow 
and identified two basics swimming states, an upstream 
oriented swinging motion around the centerline and tum- 
bling similar to passive particles [38]. In this article, we 
extend our work and analyze the motion of elongated mi- 
croswimmers, as they occur more frequently in nature, 
subject to Poiseuille flow in microchannels with different 
geometries. In our study, we will assume that gravity or 
other external forces are not present or not relevant. Our 
swimmers do not exhibit active rotation (as in Ref. |41|), 
so their local reorientation in flow is the same as for pas- 
sive elongated particles. The hydrodynamic reorientation 
in low Reynolds number flow depends on both the local 
flow vorticity and the local strain rate in the fluid. They 
both determine the periodic tumbling motion of passive 
ellipsoidal particles in flow commonly known as Jefferey 
orbit [42 . The strain rate does not influence the orienta- 
tion of spherical particles and was therefore not considered 
in our Letter [38] . 

Now, in contrast to passive particles, microswimmers 
also change their position in flow due to their active mo- 
tion. While rigid passive particles just follow streamlines 
during rotation [43j , active particles cross streamlines be- 
cause of their self-propulsion. In this article we demon- 
strate that microswimmers with elongated shape still per- 
form the characteristic swinging and tumbling motion in 
steady Poiseuille flow as observed for spherical swimmers. 
The motion is periodic in planar Poiseuille flow or in cylin- 
drical tubes with spherical cross section. Interestingly, the 
motion becomes quasiperiodic even for spherical swim- 
mers when the flow cross section becomes elliptical. 
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In the following we first introduce the equations of mo- 
tion in section |2| discuss the swimmer trajectories for dif- 
ferent flow cross sections in section [3) and conclude with 
final remarks in section [4j 



2 Equations of motion 

We study the dynamics of an elongated or prolate spheroidal 
microswimmer with aspect ratio 7 in a steady laminar 
flow Vf{r) at low Reynolds number. In bulk, the swim- 
mer moves with a constant velocity Vq = vqC^ where e 
is the swimming direction along the major axis of the 
swimmer [Fig. [lla)]. Low Reynolds number also means 
that we can neglect any inertia of the microswimmer. In 
addition, we consider the case of large rotational Peclet 
number, where rotation due to flow vorticity exceeds ro- 
tational diffusion. Furthermore, the microswimmer should 
be small compared to the lateral extension of the im- 
posed flow in the microchannel and should not disturb the 
Poiseuille flow field. Finally, we assume that the swimmer 
stays away from bounding walls so that we can neglect 
steric and hydrodynamic interactions between swimmer 
and walls which otherwise play an important role [44,45, 
I46 |l47[ l38] . Under these conditions, a passive particle sim- 
ply follows the streamlines with its center-of-mass velocity 
at position r given by v{r) = Vf{r). Its orientation moves 
on a Jefferey orbit, where the total angular velocity for a 
prolate spheroidal particle reads ^ 

flir) = ^flfir)+Gex[Eir)-e]. (1) 



(a) 



(b) 



Here we have introduced the local flow vorticity 

Hfir) = V X Vfir) (2) 
and the local strain rate 

Hr) = ^[Vvfir) + iyvfir)f]. (3) 

2_-i 

The geometry factor G = ^2^1 with G G [0, 1) depends on 
the aspect ratio 7. For spherical swimmers, where 7 = 1 
and G = 0, only the flow vorticity f2f contributes to the 
swimmer's angular velocity. An active particle with intrin- 
sic swimming speed vq and moving along its orientation e 
assumes the total velocity 



v{r) = Vf{r) + voe, 



(4) 



SO that it can easily cross streamlines. Since we do not 
assume an intrinsic rotation rate of the swimmer, the an- 
gular velocity is the same as for a passive particle [Eq. ([T])]. 
Using Eqs. ([T]) and Q the dynamics for position r{t) and 
orientation e{t) of the swimmer follows 

^r = v{r), ^e = f2{r) X e. (5) 

We consider a Poiseuille flow along the i;-direction, 

Vf{x,y) = Vf[l - (ax)^ - {byf]z, (6) 







(c) 



b=0 



® ® 

0<b<a b=a 




Fig. 1. (a) Spheroidal microswimmer with aspect ratio 7 and 
orientation e (black arrow), (b) Sketch of the Poiseuille flow 
profile in the x-z plane, (c)-(e): Sketch of the channel cross sec- 
tions (x-^-plane) with semimajor axis and semiminor axis 
and color-coded flow strength for (c) two parallel and in- 
finitely extended plates {b = 0) or a channel with (d) elliptical 
{0 < b < a) and (e) circular {b = a) cross section. 



where we use Cartesian coordinates (x, z) with the co- 
ordinate basis {x^y^z) and Vf as the maximum flow ve- 
locity in the center of the channel. For 6 = 0, Eq. ([6| 
describes the flow between two parallel and infinitely ex- 
tended plates. For b = a the flow occurs in a cylindrical 
channel, and for < 6 < a the channel has a general el- 
liptical cross section, where is the semimajor axis and 
b~^ the semiminor axis [Fig. [ijb-d)]. The flow vorticity is 
given by 

f2f{x, y) = 2vf{a^xy - b^yx) (7) 
and the strain rate is 



yyz 



yzy) 



(8) 



In the following we use rescaled units, ax ^ by ^ 
by, az and t/to t with to = {av^)'^ . We also 

introduce the dimensionless flow speed Vf = Vf/vo and 
the aspect ratio of the channel cross section, b = b/a. The 
third parameter of the system is the aspect ratio 7 of the 
swimmer. 

We note that the general configuration space for the 
swimmer dynamics is the product R^(g)S^ of position space 
(r G M^) and unit sphere (e G S^). Since we only consider 
unidirectional flow, i.e. flow in z direction, the dynamics 
of the swimmer does not depend on the z coordinate, and 
the relevant phase space becomes four-dimensional, (g) 
S^. The dimension of the relevant phase space is further 
reduced when symmetric flow geometries are used, such 
as planar or cylindrical flow. 



3 Swimmer dynamics in Poiseuille flow 

In the following we study the deterministic dynamics of 
a prolate spheroidal microswimmer in Poiseuille flow for 
several geometries. We start with the motion between two 
parallel plates {b = 0), continue with a cylindrical mi- 
crochannel with circular cross section (b — 1), and flnally 
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discuss the motion in a channel with general elliptical cross 
section {b < 1). 



3.1 Planar Poiseuille flow 



First, we discuss the motion of a microswimmer in planar 
Poiseuille flow, v f = {1 — x^)z^ with flow vorticity f = 
2vfX and strain rate E = —Vfx{xz-\-zx). For the position 
of the swimmer we employ Cartesian coordinates, r = 
XX -\- yy -\- zz^ and we expand the orientation vector e 
in the Cartesian basis using angular coordiantes ^ and O 
such that e = e^i^, 0)x + ey{0)y + e^i^, 0)z with 



cos O sin Cy = sin 0, 



- cos 6) cos IZ^. (9) 



Here ^ indicates the orientation relative to the 2:-axis and 
7^ the orientation out of the x-z-plane. We note that 
^ = 0^0 = corresponds to perfect upstream orientation 
and iZ^ = TT, = to downstream orientation, respec- 
tively. Due to translational symmetry of the flow field in 
^/-direction the relevant configuration space of the swim- 
mer reduces to three dimensions, {x, iZ^, 0} G R (g) S^. 

Using Eqs. ([5) and ([9|, the coupled equations of motion 
for (x, iZ^, 0) read 



X = — cos sin 
4^ = Vfx{l - Gcos2lZ^), 
G 

= VfX — sin 20 sin 2^. 



(10) 



Whereas the first equation describes cross-streamline mo- 
tion of the microswimmer, the second and third equations 
are the same for passive ellipsoidal particles but now the 
position X is time-dependent. When the solutions 
^{t) and 0{t) of Eqs. (10) are known, y{t) and z{t) follow 
by integration. 



y{t) = y{0) 

z{t) = z{0) 



[ sm0{t')dt\ 
Jo 

[ [Vf (1 - X{t')^) - COS 0{t') COS ^{t')] dt\ 

Jo 



3.1.1 Two-dimensional motion 

Before discussing the full three-dimensional dynamics of 
the swimmer, we consider first the two-dimensional case. 
When = 0, the swimmer only moves in the x-z-plane 
and the relevant phase space becomes two-dimensional, 
{x, ^} G M(8)S. The dynamics of position x and orientation 
^ follows from 



^ = Vfx{l - Gcos2lZ^), 
X = — sin ^. 



(12) 



For a spherical swimmer (G = 0) we recover the case 
treated in Ref. [38 . The dynamic equations are 4^ = VfX 




Fig. 2. ^-x phase portraits for Vf = 10 for a spherical swimmer 
with 7 = 1 or G = (a) and an elongated swimmer with 
7 = 5 or G = 0.923 (b). The red line is the separatrix deviding 
swinging and tumbling motion. 



and X = 
Uq sin ^ 
and cjn 



- sinlZ^", which gives the pendulum equation, ^ + 
= 0, where x plays the role of angular velocity 
— Vf. The corresponding constant of motion is the 
Hamiltonian [38] 



H -Vfx'^ + 1 - coslZ^. 
2 



(13) 



For elongated swimmers, whose dynamics is governed 
by Eqs. (12), we are also able to identify a constant of 
motion. 



Co = ^VfX^^l-go{^;G) 



(14) 



with 



^ ^ v/2G(l + G) ^ ^ 

where atanh is the inverse function of the hyberbolic tan- 
gent. We derive Cq in Appendix |A.1[ Note that Cq ^ H 
for G ^ 0. For each swimmer trajectory the constant of 
motion is uniquely defined by the initial conditions iZ^(O) 
and x{0). We note that steric interactions of the swimmer 
with the channel wall modify Gq. After contact with the 
wall at \x\ = 1 the swimmer leaves the wall with upstream 
orientation iZ^ = and the constant of motion becomes 
Go = ^ + 1 - go{0; G) ^ ^ + f (1 - |G) for small G. 

Figure [2] shows the x-^ phase portrait for ^/ = 10. 
For G = (or 7 = 1) it corresponds to the mathematical 
pendulum (a). The phase portrait for elongated swimmers 
looks qualitatively the same (b): Closed and open periodic 
trajectories exist devided by a separatrix. The trajectories 
are the curves 



x\^) = ^[go{^',G)-l^Co]. 



(16) 



The closed trajectories correspond to swinging motion of 
the microswimmer in Poiseuille flow as discussed in Ref. 
[38 . Since iZ^" = means upstream motion, the swimmer 
swings around the centerline of the channel while swim- 
ming upstream. Because of the upstream orientation, the 
flow vorticity always reorients the swimmer towards the 
centerline. After crossing the center at x = with maxi- 
mum angle iZ^max, the sign of the vorticity changes and the 
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swimmer is reoriented back towards the centerline. So, the 
swimmer experiences a periodic swinging motion around 
the centerline, in full analogy to the oscillations of a pen- 
dulum. In the channel frame, the swimmer either moves 
upstream (i < 0) for suffieciently small flow strengths Vf 
whereas it drifts downstream (i > 0) when Vf becomes 
large. Since 

z{t) = v,{t) + e,{t) = Vf[l - x{tf] - cos^{t), (17) 

the swimmer moves upstream for Vz < — or Vf < 1. 
Perfect upstream swimming in the center of the channel 
corresponds to the stable fixed point (iZ^ = 0,x = 0). 




Fig. 3. (a) The angular velocity i?^ = 1^ at position x and 
given flow strength Vf plotted versus swimmer orientation ^ 
for several shape parameters G; orange line: C = 0; light green: 
G = 0.25; dark green: G = 0.5; light blue: G = 0.75; dark blue: 
G = 1. (b): Sketch of the angular velocities of spherical and 
elongated swimmers oriented parallel and perpendicular to the 
flow. 



When the initial state of the swimmer lies outside the 
separatrix, it tumbles in the flow. The vorticity is too 
strong for the swimmer to reach the centerline. This is 
in analogy to the circling motion of a pendulum. 

Swinging and tumbling motion are devided by the sep- 
aratrix (red curves in Fig.|2|. The trajectories x{^) [Eq. ([16] 
correspond to tumbling when a position x with ^ = ir ex- 
ists or when go{7r;G) — 1 -\- Co > 0. Otherwise they are 
swinging trajectories. The limiting case is the separatrix 




Fig. 4. (a)-(b) Probability distributions P{^) for single tra- 
jectories with initial conditions ^(0) — and />(0) = 0.3 (red), 
p(0) = 0.5 (green) and p(0) = 0.8 (blue) for spherical (a) and 
elongated (b) particles with 7 = 3 (C = 0.8). (c) Angular fre- 
quencies uj for swinging and tumbling motion for t*/ = 10 for 
a spherical swimmer (orange) and an elongated swimmer with 
7 = 5 (C = 0.923) (blue) plotted versus the constant of motion 
Co. The frequency for small oscillations is cj 



G). 



which for a fixed Vf and G is given by 



x*(iZ^*; G, Vf) = ±J^ [go{^^; G) - go{7r; G)]. (18) 

V ^/ 

The separatrix possesses the constant of motion Cq = 
1 — 9q{^]G) ^ 2 — y(l — ^G). It includes the spherical 
swimmer (G = 0) with G^ = i^* = 2. For Go < G^ the 
swimmer is swinging and for Go > Cq it is tumbling. 

Linearizing the equations of motion [Eqs. ([l2|] around 
the fixed point = 0, x = 0) results in a harmonic os- 
cillator equation ^ -\- Vf{l — G)^ = with the eigen fre- 
quency ujQ = Y^^%0^^-^)- So, elongated swimmers swing 
with a smaller frequency around the centerline compared 
to spherical swimmers. The reason is that for elongated 
swimmers the reorientation rate or angular velocity is smaller 
compared to spherical swimmers, when < 7r/4. In 
particular, at fixed x the angular velocity for elongated 
swimmers [Eqs. (12), first line] depends on the orienta- 
tion, ^ = Qy = Vfx{l — Gcos2lZ>'), in contrast to spherical 
swimmers. We plot the angular velocity ^ versus ^ in 
Fig. [sja). For spherical swimmers at position x it does 
not depend on orientation ^ (orange line) while for ellip- 
soidal particles ^ varies with ^ since the strain rate of 
the flow contributes to the angular velocity. An elongated 
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Fig. 5. Swinging (a-c) and tumbling (d-e) motion of a spheroidal swimmer in flow between two parallel plates. The flow 
strength is = 10 and the aspect ratio of the swimmer is 7 = 3 (meaning G — 0.8). The initial conditions are iZ^(O) = and 
6>(0) = 7r/5 for both cases and x(0) = 0.3 for swinging and x(0) = 0.7 for tumbling motion. (a,d): Relevant three-dimensional 
phase space (x,lZ^, O). The intersection of the two constants of motion (orange) and (green) gives the periodic trajectory 
(blue intersection curve). While /^i = 2.64 is the same for swinging (a) and tumbling (d), Hi2 — 1-14 for swinging and HI2 — 2.86 
for tumbling motion. (b,e): Periodic orbits of the orientation vector e(t) which moves on the unit sphere. For swinging motion 
it is a half ellipse (b) and for tumbling motion a closed curve (e). (c,f): Swimmer trajectory r(t) in real space for swinging (c) 
and tumbling (f). The midplane is shown in grey and the trajectories are red. 



swimmer oriented parallel to the flow direction (iZ^ = or 
it) rotates with a smaller 1^ as a spherical swimmer, while 
at perpendicular orientation to the flow ^ i'7r/2) the 
angular velocity is larger. In the extreme case of an in- 
finitely thin needle {G = 1) [blue curve in Fig. |3ja)], the 
angular velocity vanishes at iZ^ = and ^ = tt and the 
needle just swims along the flow. We also sketch the re- 
orientation rate for spherical and elongated swimmers in 
Fig-ib). 

The fact that elongated swimmers rotate slower when 
oriented in flow direction has a noticeable effect on the 
orientational distribution function of the swimmer, 



P{^) (X (X (1 - Gcos2iZ^)-^(Co - 1 + ^/o(^, G))-^/2 
^ (19) 
To arrive at Eq. ( p!9|), w e have eliminated location x in 
the first line of Eqs. (|l2|) using the constant of motion Go 
of Eq. ( p!4| ). In Figs.^Ta) and (b), we show the respec- 
tive probability distributions for spherical and ellipsoidal 
swimmers. Spherical swimmers have the same probabil- 
ity distribution as the mathematical pendulum [48 . P{^) 
assumes its absolute minimum where ^ is largest. This 
occurs ai ^ = when the swimmer has its maximum dis- 
tance \x\ from the centerline [see Eqs. (12)]. In the swing- 



ing motion P{^) diverges at the centerline, where ^ = 
[green and red curve in Fig.|4ja)]. For the tumbling motion 
it stays finite since the swimmer does not reach the center- 
line [blue curve in Fig. [Zj^a)]. For ellipsoidal swimmers 4^ 
also depends on the orientation ^. Since they rotate slower 
dit ^ = when they are oriented parallel to the flow, an 
additional maximum at iZ^ = appears [Fig. |4jb)]. 



The reorientation rate ^ also determines the angular 
frequency of the swinging and tumbling motion. 



(20) 




For swinging motion, iZ^rr 
so 



^max — ^min — 



and iZ^max appear when x = 0, 



asec 



V8G/(G + l)coth ((Co - 1)V2G(1 + G) 



(21) 



where asec is the arcus secans or the inverse function of 
1/ cosx. For tumbling trajectories iZ^min = and iZVax = tt. 
We numerically integrate Eq. (20) for several geometry 
factors G and plot it versus Go in Fig. |4jc). For spher- 
ical swimmers the frequency decreases for larger oscilla- 
tions or larger Go in the swinging state in full analogy to 
the mathematical pendulum [orange curve in Fig. Qc)]. 
Larger oscillations mean that the swimmer orients more 
towards a perpendicular orientation relative to the flow 
direction. Since for elongated swimmers the reorientation 
rate is largest at = 7r/2, the frequency uj first increases 
until it reaches a maximum for a specific swinging state de- 
termined by a specific value of Go [blue curve in Fig.Qc)]. 
In the tumbling state uj increases with Go for both the 
spherical and elongated swimmer. 
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Fig. 6. Swinging and tumbling motion for a microswimmer 
in a cylindrical tube for the flow strength Vf — 10 and swim- 
mer aspect ratio 7 = 2 (C = 0.6). The initial conditions are 
iZ/(0) = 7r/4, p(0) = 0.3 and O{0) = 7r/5 for swinging mo- 
tion and ^{0) = 0, p(0) = 0.7 and O{0) = 7r/5 for tumbling, 
(a) Typical trajectories for swinging (orange) and tumbling 
(black) motion in the p-^-O phase space. The trajectories r(t) 
in real space are shown in (b) for tumbling motion and in (c) 
for swinging motion together with projections onto the channel 
cross section. Note that the z direction is squeezed relative to 
the lateral dimension. 



3.1.2 Three-dimensional motion 



For the three-dimensional motion [Eqs. (10)] we a re ab le 
to identify two constants of motion (see Appendix A.2), 



Ki =tan^ 0/(1 -G cos 2lZ^), 
K2 = \vfx^ -gi{^,0',G) 



with 



<7i(!^,0;G) 



at anh ( ^J2G/{l + G) coslp-cos©) 
\/2G(l + G) • 



(22) 



(23) 



The intersection of both constants and in the x- 
^-0 phase space defines again periodic solutions of the 
ellipsoidal swimmer in the planar Poiseuille flow. For = 
0, hzi vanishes and gives Cq. We thus recover the two- 
dimensional motion discussed in the previous section. In 
the full three-dimensional case, we can again divide the 
periodic solutions into swinging and tumbling motions. 

Figure [5] illustrates specific examples for Vf = 10 and 
G = 0.8. Figure |5|a) shows the phases space curve for a 
typical swinging motion and Fig. jsjd) for the swimmer 
in a tumbling state. The initial conditions are ^{0) = 
and 0(0) = 7r/5 for both cases, and x(0) = 0.3 for (a) 
and x{0) = 0.7 for (d). So the constant of motion hzi is the 
same for both cases and it is illustrated by the orange sur- 
faces. The second constant of motion (green surfaces) 
depends on x and thereby defines if swinging or tumbling 
occurs. The intersection between a^i and hi2 is the actual 
trajectory in phase space (blue curve). Figures jsjb) and 
(e) show the motion of the orientation vector e(T). Since 
it is a unit vector, its arrowhead moves on the unit sphere 
S^. For the swinging motion e moves back and forth with 



e^y > on a trajectory that resembles a half ellipse (b), 
while it circles around a closed trajectory during tumbling 
(d). The closed trajectory is similar to the Jefferey orbit of 
tumbling passive particles 142] , however for active particles 
X is not constant. In Fig. pjc) and (f) the trajecories r{t) 
are shown. For swinging motion the middle plane {x = 0) 
is crossed periodically (c), while tumbling motion always 
occurs on either side of the middle plane (f). Both trajec- 
tories look two-dimensional tilted against the z direction. 
Since the component ey{t) oscillates with time, they are 
clearly three-dimensional. 

However, for spherical swimmers G = and = 
[Eq. (10)], so that the component Cy = sin0 is constant in 
time. The equation of motion for ^ is again the Pendulum 
equation but it depends on the initial orientation Cy = 
const., 

^ ^^Vfsin^ = 0, (24) 



where ^ = y 1 — = const is the length of the orien- 
tation vector in the x-z plane. A swimmer starting with 
ey ^ swings or tumbles at smaller frequency ujq = y^^^ 
compared to a swimmer with = on a two-dimensional 
trajectory tilted against the z direction. So, the position 
y{t) increases linearly in positive direction for 0(0) > 
or in negative direction for 0(0) < 0. Since = const 
the relevant phase space is the x-^-spdiCe similar to the 
two-dimensional problem. 




Fig. 7. (a)-(c) Two- frequency quasiperiodic swinging motion 
for the initial conditions ^(0) = 0, cp{0) = 0, p(0) = 0.3, 
0(0) 7r/5 and parameters 6 = 2, 7 = 2, = 10. (d)-(f) 
Tumbling motion for the initial conditions ^(0) — 0, (p{0) — 
0, p{0) 0.7, 0(0) 7r/5 and parameters b 1.4, 7 = 2 
and Vf = 10. (a) and (d): Three-dimensional motion in the 
elliptical channel, (b) and (e): Trajectories projected onto the 
cross section of the channel, (c) and (f): Time evolution of the 
swimmer orientation component Czit). 



3.2 Swimming in a cylindrical tube 

Now we discuss the motion of a swimmer in a cylindrical 
microchannel. We use cylindrical coordinates p^(p,z and 
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(a) (b) (c) (d) 




Fig. 8. Projected phase space and Poincare sections for quasiperiodic swinging (a)-(d) and tumbling (e)-(li) motion. The initial 
conditions and system parameters are the same as in Fig. [7] The projections of the trajectories onto the p-^-O phase space lie 
on two-dimensional surfaces for swinging (a) and tumbling (e) motion, (b) and (f): same trajectories as in (a) and (e) but with 
the coordinate cp color-coded, (c) and (g): Poincare sections at cp = (cp > 0). One-frequency behavior in the Poincare sections 
indicate two- frequency motion in phase space ^49j. (d) and (h): Projections onto the p-^ phase space where cp is color-coded. 



express the orientation vector e in the cylindrical coor- 



dinate basis {p^(p^z)^ 
ez{^,9)z with 



Cp = — cos O sin iZ^, = sin 0, = — cos cos ^. (25) 

Because of translational symmetry in z direction and ro- 
tational symmetry with respect to (p, only the equations 
for p, ^ and are coupled, 

p = — cos sin iZ^, 

4^ = Vfp[l -Gcos2^] - sin0tan0coslZ^/p, (26) 
(9 = (l + 2Gvfp^ cos cos ^) sin sin ^/p. 

The relevant phase space R(8)S^ is again three-dimensional. 
For the spherical microswimmer (G = 0) we already iden- 
tified two constants of motion. They lead to periodic tra- 
jectories with helical or helical-like swinging and tumbling 
motion [38^. We derive the constants of motion in Ap- 
pendix |A.3 We were not able to construct constants of 
motion for an elongated swimmer in cylindrical Poiseuille 
flow. Nevertheless, we numerically find that it performs 
again periodic motions similar to the spherical swimmer. 
Figure [6]^a) shows two closed phase space orbits, both for 
the swinging and tumbling state, for a swimmer with as- 
pect ratio 7 = 2 or geometry factor G = 0.6. We also il- 
lustrate the three-dimensional trajecories in the tube and 
their two-dimensional projections onto the channel cross 
section Figs. l6Fb) - (c). For the swinging motion we have 



chosen initial conditions such that the projection onto the 
cross-sectional plane is periodic. 



3.3 Quasi-periodic motion in a channel with elliptical 
cross section 

Now we consider the motion in a channel with elliptical 
cross section (h > 1). We again use a cylindrical coordinate 
system as for the cylindrical channel. The symmetry in (p 
direction no longer exists and the relevant phase space is 
four-dimensional, (p, iZ^, 0, (p) G S^. The equations of 
motion for the relevant coordinates read 



p = — cos sin 
4^ = Vfp{l - Gcos2iZ^) 
— sin tan cos ^/ p 



■ cos 2(p 







Vfp 



sin sin ^ 



VfP 



(1 + G) sin 2(p sin iZ^ tan 0, 



[l + 2Gvfp^ cos cos lZ>'(cos^ + 6^ sin^ Lp) 



(l-G cos 20) cos sin 2ip, 



(p = sm0/ p. 



(27) 
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We again solve them numerically and determine the swim- 
mer trajectories. Interestingly, both spherical and elon- 
gated swimmers now follow quasiperiodic trajectories. Fig- 
ure [7] shows trajectories for a quasiperiodic swinging [Fig. 
[rj^aj] and tumbling motion [Fig. ^d)] and their projec- 
tions onto the elliptical channel cross sections [Figs, 
and (e)]. We also show the time evolution of the z com- 
ponent of the orientation vector, ez{t)^ which shows two- 
frequency quasiperiodic behavior [Figs.[7|^c) and (f)] as we 
will also demonstrate in Fig.jSj This indicates that the dy- 
namics of a single trajectory takes place on a torus in the 
four-dimensional phase space. In Fig. [S] we further illus- 
trate the complex dynamics of the quasiperiodic swing- 
ing [Figs. (Sja) - (d)] and tumbling [Figs, (sje) - (h)] dy- 
namics for the same parameters as in Fig. [7[ Since the 
quasiperiodic trajectories lie on a torus, they define a 
two-dimensional surface when projected onto the three- 
dimensional subspace of the full phase space. Fig- 
ures (sja) and (e) illustrates these complex surfaces. In 
addition, we plot their Poincare sections for (p = 0, cp > 
in Fig. Isjc) and (g). The resulting closed curves confirm 
the two- frequency quasiperiodic motion on a torus [49 . 
In Fig. [sjb) and (f) we show the surfaces again with the 
azimuthal angle cp color-coded. For the swinging motion 
[Fig. [sjb)] one can clearly identify specific regions in the 
space for each value of (p. The projection of the 
surface onto the p-^ plane [see Fig. [sj^d)] resembles the 
closed trajectories of the swinging motion discussed ear- 
lier but now the curve is smeared out in the p-^ plane due 
to the quasiperiodicity. The same is true for the tumbling 
motion, which we illustrate further in Figs, [sj^f) and (h). 
The projected phase space trajectory in the p-^ plane is 
a smeared-out version of previous examples (Figs. [sj^h)). 
However, the azimuthal angle Lp in the color-coded surface 
of Figs.jsjf) is distributed over the whole phase-space tra- 
jectory. 

Finally we note that we also studied the dynamics of 
spherical swimmers which also turned out to be quasiperi- 
odic. 



4 Conclusion 

To conclude, we studied the two- and three-dimensional 
deterministic dynamics of elongated swimmers in Poiseuille 
fiow at low Reynolds number generalizing our work from 
Ref. [38]. The state of the system is given by five vari- 
ables, three for the position and two for the orientation 
in the channel, respectively. For all swimmers we find the 
same characteristic swimming states. Depending on initial 
conditions and parameters - aspect ratio of the swimmer 
(7), geometry of the channel cross section (6), and fiow 
strength {vf) - the microswimmer either swings with up- 
stream orientation around the centerline of the channel 
or it tumbles at sufficiently large fiow vorticities. The fre- 
quency of swinging and tumbling motion depends on the 
aspect ratio, in particular, it is smaller for elongated com- 
pared to spherical swimmers. For the motion between two 
infinitely extended parallel plates and for the motion in a 



cylindrical channel the relevant phase space is only three- 
dimensional due to the symmetry of the channel cross sec- 
tion. In these cases we were able to identify two constants 
of motion which results in a periodic motion visible in the 
three-dimensional phase space. Interestingly, when we re- 
duce the symmetry by choosing a channel with elliptical 
cross section, the phase space becomes four-dimensional 
and quasiperiodic motion occurs. 

Note that in the framework of dynamical systems the 
motion of swimmers in Poiseuille fiow is not considered 
as dissipative. In laminar fiow noninter acting swimmers 
move on separate trajectories without aggregating, their 
dynamics does not reach a common fixpoint or limit cycle. 
This changes when hydrodynamic interactions between a 
swimmer and the channel wall play a significant role as 
demonstrated in Ref. [38 where we explicitely identified 
fixpoints and limit cycles for the spherical swimmer. Other 
examples are external forces such as gravitation j34H35( l5Q] 
or external light sources [51 that infiuence the reorienta- 
tion of the microswimmer s. 



A Constants of motion 
A.l 2D planar flow 

A constant of motion can be found by deviding the first 
equation by the second equation of Eqs. (12), 



X 



dx 



sinlZ^" 



This expression can be separated and integrated, 

sin ^d"^ 



G cos 2^ ' 



leading to 



2 



atanh (^^2G/{l + G) cos^^ 



(28) 



(29) 



(30) 



where c is a constant and Co 
motion defined in Eq. (14). 



V2G(1 + G) 

c + 1 is the constant of 



A. 2 3D planar flow 

The first constant of motion i^i is identified by calculating 
d0 Gsm20sm2^ 



4r 2 1 - Gcos2iZ^* 

We again separate the variables and integrate, 
d0 r sm2^d^ 



Gs'm20 / l-Gcos2!Z^' 



(31) 



(32) 



resulting in 
1 



2^[log(l-Gcos2lZ^)]+c= - [log(tan0)]. (33) 
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Taking the exponential of this expression and defining ki 
as Hii = e~^^ results in the first constant of motion in 



Eqs. (22). Using /^i, eliminating and calculating 

- = ^ = (1 - Gcos2iZ^) JlT^^^^l^^G^^^, 

(34) 

after separation of variables and integration then leads to 

(35) 



Vf 2 _ atanh(A/2G/(l + G) cos ^ cos 0) 



V2G(1 + G) 



where we defined as = 1 + c as the second constant 
of motion in Eqs. (22). 



A.3 Cylindrical flow for G = 

The equations of motion for spherical particles swimming 
in a cylindrical Poiseuille flow are given by Eqs. (26) but 



with G = 0. Using Eqs. ( |26| ) the first constant of motion 
for a spherical swimmer is determined by taking 



= — tan6)/p, 



(36) 



o_de 

p dp 

after separation of variables and integration leading to 

- log(sin 0) = log(p) + c, (37) 

and to the first constant of motion L = = psinO 
which is proportional to the angular momentum of the 
swimmer in the z direction. The second constant of motion 
is identified by taking 



dO 



tan0 



sin'^ O sin ^ tan ^ 



(38) 



This nonlinear differential equation can be transformed 
into a linear equation by using the ansatz w = cos ^ lead- 
ing to 

dw 

de 

The solution is 



; tan O 



(39) 



w 



COS ^ 



2 sin O tan O cos O 
and we define the second constant of motion as 



(40) 



M 



1 



1 



L^Vf + 1 = + 1 - cos cos ^. (41) 
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